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Abstract. In modeling physical systems, it is sometimes useful to construct 
border bases of O-dimensional polynomial ideals which are contained in the 
ideal generated by a given set of polynomials. We define and construct such 
subideal border bases, provide some basic properties and generalize a suitable 
variant of the Buchberger-Moller algorithm as well as the AVI-algorithm of [5] 
to the subideal setting. The subideal version of the AVI-algorithm is then 
applied to an actual industrial problem. 
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1. Introduction 

In [5] an algorithm was introduced which computes an approximate border basis 
consisting of unitary polynomials that vanish approximately at a given set of points. 
It has been shown that this AVI-algorithm is useful for modeling physical systems 
based on a set of measured data points. More precisely, given a finite point set 
X = {pi, , . . ,p s } C [—1,1]™, the AVI-algorithm computes an order ideal O of 
terms in T™ and an O -border prebasis G = {g\, . . . , g v } such that 

(1) the unitary polynomials <7j/||<?i|| vanish e -approximately at X, where e > 
is a given threshold number, and 

(2) the normal remainders of the S-polynomials S(gi,gj) for gi,gj with neigh- 
boring border terms are smaller than e . 

Abstractly speaking, the last condition means that the point in the moduli space 
corresponding to G is "close" to the border basis scheme (see [IT] and [8]). In 
practical applications, the AVI-algorithm turns out to be very stable and useful. 
With a judicial choice of the threshold number e, it is able to discover simple 
polynomial relations which exist in the data with high reliability. For instance, 
it discovers simple physical laws inherent in measured data without the need of 
imposing model equations. 
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However, in some situations physical information may be available which is not 
contained in the data points X or we may have exact physical knowledge which is 
only approximately represented by the data points. An example for this phenom- 
enon will be discussed in Section 6. For instance, we may want to impose certain 
vanishing conditions on the model equations we are constructing. Using Hilbert's 
Nullstellensatz this translates to saying that what we are looking for is the inter- 
section of the vanishing ideal of X with a given ideal J C R[a;i, . . . , x n ] whose 
generators represent the vanishing conditions we want to impose. 

In order to be able to deal with this approximate situation, it is first necessary to 
generalize the exact version of the computation of vanishing ideals to the subideal 
setting. Then this theory will serve as a guide and a motivation for the approxi- 
mate case. Therefore this paper begins in Section 2 with the definition and basic 
properties of subideal border bases. 

Given a O-dimensional ideal 7 in a polynomial ring P = K[x±, . . . , x n ] over a field 
and a set of polynomials F = {/i, . . . , f m } generating an ideal J = (F) , a subideal 
border basis of I corresponds to a set of polynomials Of = 0\ ■ /i U • • • U O m ■ f m , 
where the Oi are order ideals of terms, such that the residue classes of the elements 
oi Op form a K -basis of J /(ID J) = (I + J)/I . Clearly, this generalizes the case 
F = {1}, i.e. the "usual" border basis theory. We show that subideal border bases 
always exist and explain a method to construct them from a border basis of /. 
Moreover, we discuss some uniqueness properties of subideal border bases. 

The foundation of any further development of the theory of subideal border bases 
is a generalization of the Border Division Algorithm (see [10], 6.4.11) to the subideal 
case. This foundation is laid in Section 3 where we also study higher Of -borders, 
the Of -index, and show that a subideal border basis of I generates I fl J . 

In Section 4, we generalize the Buchberger-Moller algorithm (BM-algorithm) for 
computing vanishing ideals of point sets to the subideal setting. More precisely, 
we generalize a version of the BM-algorithm which proceeds blockwise degree-by- 
degree and produces a border basis of the vanishing ideal. Similarly, the subideal 
version of the BM-algorithm (cf. Algorithm 14. 2p computes an O a (Ix)f -subideal 
border basis of 1%, where O a (I\) is the complement of a leading term ideal of the 
vanishing ideal Ix of X. 

Next, in Section 5, we turn to the setting of Approximate Computational Alge- 
bra. We work in the polynomial ring . . . , x n ] over the reals and assume that 
X C [—1, 1]™ is a finite set of (measured, imprecise) points. We define approximate 
Of -subideal border bases and generalize the AVI-algorithm from [5], Thm. 3.3 to 
the subideal case. 

Let us point out that the subideal version of the AVI algorithm contains a sub- 
stantial difference to the traditional way of computing approximate vanishing ideals, 
e.g. as in [I]. Namely, the AVI algorithm produces a set of polynomials which van- 
ish approximately at the given data points, but we do not demand that there exists 
a "nearby" set of points at which these polynomials vanish exactly. The latter 
requirement has turned out to be too restrictive for real-world applications, for 
instance the one we explain in the last section. There we provide an example for 
the application of these techniques to the problem of production allocation in the 
oil industry. 

Unless explicitly stated otherwise, we use the notation and definitions of [9] 
and [10]. We shall assume that the reader has some familiarity with the theory of 
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exact and approximate border bases (see for instance [5], [6], [7], [8], Section 6.4 
of [10], and [12]). 

2. Subideal Border Bases 

Here we are interested in a "relative" version of the notion of border bases in 
the following sense. Let if be a field, let P = K\x\, . . . , x n ] be a polynomial ring, 
let T" be its monoid of terms, let O be an order ideal in T n , and let I C P be a 
O-dimensional ideal. 

Suppose we are given a further polynomial ideal J = (/i, . . . , f m ) of P , where 
F = {/i, . . . , / m } C P \ {0} . Our goal is to describe and compute the intersection 
ideal I D J as a subideal of J. By Noether's isomorphism theorem, we have 
J I (I n J) = (I + J) 1 1 C P/I ■ Therefore J has a finite if -vector space basis 
modulo I n J . Now we are looking for the following special kind of vector space 
basis. 

Definition 2.1. Let O be an order ideal of terms in T™ whose residue classes form 
a if -vector space basis of P/I . 

(1) For i — l,...,m, let (9, C O be an order ideal. Then the set Of = 
01 • /l U • • • U O m ■ f m is called an F -order ideal. Its elements, i.e. products 
of the form tfi with t d Oi will be called F -terms. 

(2) If Of = 0\ ■ fi U • ■ ■ U O m ■ f m is an i^-order ideal whose residue classes 
form a A" -vector space basis of J /(I flJ), we say that the ideal I has an 
Of -subideal border basis. 

Notice that an F -term may be viewed as a generalization of the usual notion 
of term by using F = {1}. Similarly, F -order ideals generalize the usual order 
ideals. It is natural to ask whether every ideal I supporting an O -border bases has 
an Op -subideal border basis for some F -order ideal Of- The next proposition 
answers this positively. 

Proposition 2.2. Let I C P be a O-dimensional ideal, and let J = (/i, . . . , f m ) C 
P be any ideal. 

(1) Given an order ideal OcT™ whose residue classes generate the K -vector 
space P/I , there exists an order ideal O C O whose residue classes from 
a K -vector space basis of P/I . 

(2) Let U C P m be a P -submodule, and let Ox,..., O m be order ideals in T™ 
such that the residue classes of 0\ex U ■ ■ • U O m e rn generate the K -vector 
space P m /U . Then there exist order ideals Oi C Oi such that the residue 
classes of Oxex U • • • U O m e m form a K -vector space basis of P m /U . 

(3) If O is an order ideal whose residue classes from a K -vector space basis 
of P/I then there exist order ideals Oi C O such that the residue classes 
of Of = Ox fx U • • • U O m f m are a K -vector space basis of J /{I O J) . In 
other words, the ideal I has an Of -subideal border basis. 

Proof. First we show (1). We construct the order ideal O inductively. To this end, 
we choose a degree compatible term ordering a and order O = {tx, • ■ ■ , t^} such 
that tx <cr ' ' ' <cr . In particular, we have tx = 1. Since I C P , we can start by 
putting tx into O and removing it from O . 
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For the induction step, we consider the er-smallest term t% which is still in O . 
If the residue class of U in P/I is if -linearly dependent on the residue classes of 
the elements in O , we cancel U and all of its multiples in O . Each of these terms 
can be rewritten modulo J as a linear combination of smaller terms w.r.t. a. If 
the residue class of ti in P/I is if -linearly independent of the residue classes of 
the elements in O, we append t{ to O and remove it from O . In this way, the 
residue classes of the elements of O are always if -linearly independent in P/I, 
and every element of O can be rewritten modulo I as a if -linear combination of 
the elements of the final set O . 

Now we show (2). Let U be the idealization of U in P = if [x\, . . . , x n , ei, . . . , e m ] 
(see [10], Section 4.7.B). The residue classes of the elements of the order ideal 
O = Oi ei U •• • U O m e m generate the if -vector space P/U . Now it suffices to 
apply (1) and to note that every order subideal of O has the indicated form. 

Finally, we prove (3). Consider the P-linear map ip : P m — > P/I defined by 
ch > f% + I ■ Its image is the ideal (i + J)/I . Let U = kei-(p) . Then tp induces an 
isomorphism (p : P m /U = (7 + J)/ 1. To get a system of generators of P m /U , it 
suffices to find a system of generators of the ideal (i + J) /I = (fi+ 1, ■ . ■ , f m + 1) ■ 
As a vector space, this ideal is generated by O ■ (/i + I) U • • • U O ■ (f m + I) . 
The preimages of these generators are the elements of Oe\ U • ■ ■ U Oe m . Now an 
application of (2) finishes the proof. □ 

Based on this proposition, we can construct an F-order ideal such that a given 
ideal I has an Of -subideal border basis. The following example illustrates the 
method. 

Example 2.3. Let P = Q[x,y], let I = (x 2 — x, y 2 — y) , let O = {1, x, y, xy} , and 
let J = (x + y) . Then I has an O -border basis and therefore also an Op -subideal 
border basis w.r.t. F = {/} for / = x + y . 

To construct a suitable F -order ideal, we start with Op = {1 • /}. Then we 
put x ■ f and y ■ f into O p , since we have x ■ f = xy + x and y ■ f = xy + y 
modulo i, and since {x + y, xy + x, xy + y} is Q-linearly independent in P/I . 
Next xy ■ f = 2xy = x ■ f + y ■ f — 1 • / implies that we are done. The result is 
that Of = {/, xf, yf} is an F -order ideal for which I has an Oi?-subideal border 
basis. 

At this point it is time to explain the choice of the term "subideal border basis" 
in the above definition. 

Definition 2.4. Let F = {/i, . . . ,/ m } C P\{0}, and let O f = 1 f 1 U---UO m f m 
be an F-order ideal. We write Of = {tif ai , ■ ■ ■ ,t^,f a } with an € {1, . . . ,m} and 

U eO«,. 

(1) The set of polynomials OOf — {x\O f U • • • U \ Of is called the 
border of Of ■ 

(2) Let OOf = {bifp i: ■ ■ ■ ,b v fp„}. A set of polynomials G = {gi, . . . ,g„} is 
called an Of -subideal border prebasis if gj = bjfp j — Y]f—i Ci^ifcm with 
cy, . . . , c p j e if for j = 1, . . . , v . 

(3) An Of -subideal border prebasis G is called an Of -subideal border basis 
of an ideal i if G is contained in i and the residue classes of the elements 
of Of form a if -vector space basis of J /(I n J) . 
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In this terminology, the last part of the preceding proposition can be rephrased 
as follows. 

Corollary 2.5. Let O be an order ideal in T™, and let I C P be a 0- dimensional 
ideal which has an O -border basis. Then I has an Of -subideal border basis for 
every ideal J = (fx,..., f m ) and F = {ft, f m } C P \ {0} . 

In the setting of Example 12.31 the Of -subideal border basis of / can be con- 
structed as follows. 

Example 2.6. The border of the .F-order ideal Of = {/, xf, yf} is OOf — 

{x f,xyf,y 2 f}. We compute modulo / and find x 2 f = xf , xyf = xf + yf — f , 
and y 2 f = yf . Therefore the set G = {x 2 f - xf, xyf - xf -yf + f, y 2 f - yf} is 
an Of -subideal border basis of /. 

If an ideal has an Of -subideal border basis, the elements of this basis are 
uniquely determined. This follows exactly as in the case J = (1) , i.e. the case 
of the usual border bases (see [TU], 6.4.17 and 6.4.18). Notice, however, that a set 
of polynomials may be an F-order ideal in several different ways. This is illustrated 
by the following example. 

Example 2.7. Let P = Q[x,y], let / = (x 2 — x, y 2 — y) , and let J = (x, y) . 
Clearly, the ideal / has an -border basis for O = {1, x, y, xy} , namely the set 
G = {x 2 — x, x 2 y — xy, xy 2 — xy, y 2 — y\. Hence the ideal / also has an Of -subideal 
border basis for F — {x, y} . Here we can use both Of — {1, y} ■ x U {1} • y and 
F = {l}-xU{l,x}-y. 

This example shows also another phenomenon: an F-term can simultaneously 
be contained in Of and in dOp ■ For instance, if we use Of = {l,y}-xi){l}-y , the 
term xy is both contained in {l,y} ■ x and in the border of {1} • y. The resulting 
subideal border basis will contain the polynomial xy — xy = . 

Finally, we give an example where a term is in OOf in two different ways, so 
that a subideal border basis polynomial is repeated. 

Example 2.8. Let I = (x 2 — x, y 2 — y, xy) C Q[x,y], and let J = (x, y) cQ[i,j]. 
Then the subideal border basis of / with respect to Of — {1} • x U {1} • y is 
G = {x 2 — x, xy, xy, y 2 — y} where xy appears both in <9{1} • x and in d{l} ■ y . 

3. The Subideal Border Division Algorithm 

A central result in the construction of any Grobner-basis-like theory is a suit- 
able version of the division algorithm (for the classical case, see for instance [5], 
Thm. 1.6.4 and specifically for border bases, see [10], Prop. 6.4.11). Before we can 
present a subideal border basis version, we need a few additional definitions. 

Definition 3.1. Let F = {fx, f m } C P\ {0} , and let O f be an F-order ideal. 

(1) The first border closure of O f is dO^ = O f U dO F ■ 

(2) For every k > 1, we inductively define the (k + l) st border of Of by 
d k+1 F = d{d k O F ) and the (k + l) st border closure of O f by d k+1 F = 
d k O F Ud k+1 F . 

(3) Finally, we let d O F = d O F = O f ■ 

Using these higher borders, the set T"/iU---UF f m is partitioned as follows. 
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Proposition 3.2. Let F = {/i,...,/ m } C P\ {0} , and let Of be an F -order 
ideal. 

(1) For every k > 0, we have a disjoint union 

(2) For every k > 0, we have d k Op — T" fc ■ Of ■ 

(3) For every k > 1 , we have d k O F = T£ • O f \ T n <k ■ O f ■ 

(4) We have (J i=0 T n • fi = lj°^ d 3 Of , where the right-hand side is a disjoint 
union. 

(5) Any F -term tfi G T n • fi \ Op is divisible by an F -term in dOp ■ 

Proof. First we show (1) by induction on k. For k = 0, the claim follows from 
the definition. For k = 1 , we have d 1 Q F = d a O F U d 1 F by Definition O a. 
Inductively, it follows that d k+1 G F = d k O F U k+1 O F = Ui=o dl °F ■ Thi s is a 
disjoint union, since d k+1 F n d k O F = in each step. 

Next we prove claim (2). Again we proceed by induction on k, the case k = 
being obviously true. Inductively, we have k+1 O F = d k O F Ud k+1 F = T" fc • 
O f U T™ • (T^ fe • Op) = Tl k+1 ■ Op . 

Claim (3) is a consequence of (2) and the equality d k Op = d k Op\d k ~ 1 Op . The 
fourth claim follows from the observation that, by (2), every F'-term is in d k Op 
for some k > 0. 

Finally, claim (5) holds because (4) implies that tfi G d k OF for some fe > 1, 
and by (3) this is equivalent to the existence of a factorization t = t't" where 
deg(i') = k - 1 and t" fi G dO F . □ 

In view of this result, the following definition appears natural. 

Definition 3.3. Let F = {fx, f m } C P\ {0} , and let Of be an F-order ideal. 

(1) For an F-term tfi G Op , we define ind 0j? (t/;) = min{fc > | tfi G d k O F } 
and call it the Op -index of tfi. 

(2) Given a non-zero polynomial / G J, we write / = pifi + ■ ■ ■ +p m fm with 
Pi G P and we let V = (pifi, . . . ,p m f m ) ■ Then the number 

indo F (V) = max{indo F (i/j) | i G {1, . . . , m}, i G Supp(pi)} 

is called the O f -index of the representation V of /. 

In other words, the Op -index of t fi is the unique number k > such that 
i/j G d k O F . Note that the O^-index of a polynomial / G J depends on the 
representation of / in terms of the generators of J. It is not clear how to find a 
representation V which yields the smallest inde^T 3 ). Using the Subideal Border 
Division Algorithm, we shall address this point below. 

The following proposition collects some basic properties of the Op -index. 

Proposition 3.4. Let F — {/i, . . . , f m } C P \ {0}, and let Op be an F -order 
ideal. 

(1) For an F-term tfi G T" • fi, the number k = inde> F (i/j) is the smallest 
natural number such that there exists a factorization t — t't" with a term 
t' G T" of degree k and with f'fceOp. 

(2) Given t G T" and an F -term t' fi G T" • fi , we have 

\n&o P {U'fi) < deg(t)+ind OF (t'fi). 
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(3) For /, g € J\{0} such that f + g ^ , we write f = pi/i + ••• +p m fm and 

g = qih H 1" <?m/m wift Pi, ?j € P, and we Zei "P = (pi/i, . . . ,p m fm) 

and Q = (qifi, . ■ ■ , q m fm) ■ Then we have 

ind OF (P + Q) < max{ind OF (P), ind GF (Q)}. 

(4) Given f E J \ {0} , me write / = pi/i + • ■ ■ + p m fm with pi E P and Zei 

= (Pi/i) ■ ■ ■ ,Pmfm) ■ For every g E P \ {0} , we f/iera ftaae 

in& OF (gV) < deg(g) + md 0f ,(P), 

Proof. The first claim follows from Prop. 13.21 The second claim follows from the 
first. The third claim is a consequence of the fact that every P-term appearing in 
V + Q appears in P or Q. The last claim follows from (2) and the observation 
that gV is a if -linear combination of tuples tV with t E Supp(o) . □ 

Now we have collected enough material to formulate and prove the subideal 
version of the Border Division Algorithm. 

Algorithm 3.5. (The Subideal Border Division Algorithm) 

Let F = {fi,...,f m } C P\{0} 7 let O f = {t 1 f ai ,...,t fl f a J be an F -order 
ideal where ct; E {1, . . . , m} and ti E O ai , let dOp = {b\fp 1 , . . . , b v fpj\ be its 
border, and let {gi, . . . ,g v } be an Op -subideal border prebasis, where gj = bjfp. — 
Yji=i °ijtifoti with c\j, . . . , Cpj E K for j = 1, . . . , v . Given a polynomial f E J , 
we write f = pi/i + • • • + p m f m and consider the following instructions. 

Dl Let h\ = ■ ■ ■ = h v = 0, c\ = ■ ■ ■ = c M = 0, and Q = {qxfx, . . . , q m f m ) with 

qi=Pi for i= 1, ...,m. 
D2 // Q = (0, . . . , 0) then return (hi, . . . , h v , C\, . . . , c„) and stop. 
D3 // inde> F (<2) = then find Ci, . . . , c M E K such that qifi + ■ • ■ + q m f m = 

cihf ai H h c^t^fa^ . Return (hi, . . . , h u , c x , . . . , c M ) and stop. 

D4 // indo F (Q) > then determine the smallest index i E {l,...,m} such 
that there exists a term t E Supp(gi) with indo F (tfi) — inde> F (Q) . Choose 
such a term t . Let a E K be the coefficient of t in qi . Next, determine the 
smallest index j G {1, . . . , v\ such that t factors as t — t' t" with a term t' 
of degree mdo F (tfi) — 1 and with t" fi — bjfp, E dOp . Subtract the tuple 
corresponding to the representation 

a t' gj — at' bjfp 3 - °ij a hfa, 

i-l 

from Q, add at' to hj , and continue with step D2. 

This is an algorithm which returns a tuple (hi, . . . , h v , ci, . . . , c M ) E P v x such 
that 

f = higi H h h u g u + cihf ai H h c^t^fa^ 

and deg(hi) < ind OF ("P) - 1 for V = (pi/i, • • . ,p m fm) an d for all i € {1, . . . , v} 
with hi 7^ . This representation does not depend on the choice of the term t in 
step D4. 

Proof. First we show that all steps can be executed. In step D3, the condition 
indo F (Q) = implies that all F -terms tfi with t E Supp(oi) are contained in O f . 
In step D4, the definition of indo F (Q) implies that a term t of the desired kind 
exists. By Proposition 13.41 1. this term t has a factorization t = t't" with the 
desired properties. 
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Next we prove termination by showing that step D4 is performed only finitely 
many times. Let us investigate the subtraction of the representation of at' 'gj 
from Q. By the choice of t' , the OF-hidex of t'bjfp. is deg(i') more than the 
Of -index of bjfp. . By Prop. [3~4l b. this is the maximal increase, and the Of -index 
of the other F -terms in the representation of at' gj is smaller than indo jr (Q) . Thus 
the number of F -terms in Q of maximal O^-hidex decreases by the subtraction, 
and after finitely many steps the algorithm reaches step D2 or D3 and stops. 

Finally, we prove correctness. To do so, we show that the equality 

/ = <7i/i H 1- <7m./m + h x g x -\ V h v g v + Citi/ Ql H h c^/^ 

is an invariant of the algorithm. It is satisfied at the end of step Dl. The constants 
ci, . . . , c M are only changed in step D3. In this case the contribution qifi + ■ ■ ■ + 
qmfm t° the above equality is replaced by the equal contribution Citif ai + • ■ • + 
Cfitfifa^ ■ The tuple Q is only changed in step D4. There the subtraction of the 
representation of at'gj from Q and the corresponding change in qifi + ■ ■ ■ + q m f m 
are compensated by the addition of at' to hj and the corresponding change in 
higi + ■ • ■ + h v g v . When the algorithm stops, we have q± = ■ ■ ■ = q m = Q . This 
proves the claimed representation of /. Moreover, only terms of degree deg(i') < 
inde> F (Q) — 1 < indo F ('P) — 1 are added to hj . 

The additional claim that the result of the algorithm does not depend on the 
choice of t in step D4 follows from the observation that tfi is replaced by F -terms 
of strictly smaller O p -index. Thus the different executions of step D4 correspond- 
ing to the reduction of several F -terms of maximal Of -index in Q do not interfere 
with one another, and the final result - after all those F-terms have been rewritten 
- is independent of the order in which they are taken care of. □ 

Notice that in step D4 the algorithm uses a term t which is not uniquely de- 
termined. Also there may be several factorizations of t. We choose the indices i 
and j minimally to determine this step of the algorithm uniquely, but this particu- 
lar choice is not forced upon us. Moreover, it is clear that the result of the division 
depends on the numbering of the elements of 80f- 

As indicated above, the Subideal Border Division Algorithm has important im- 
plications. The following corollaries comprise a few of them. 

Corollary 3.6. (Subideal Border Bases and Special Generation) 

In the setting of the algorithm, let I = (G) . Then the set G is an Of -subideal 
border basis of I if and only if one of the following equivalent conditions is satisfied. 

(Ai) For every non-zero polynomial f G ID J with a representation f — pxfi + 
■ • -+Pmfm and V = (pi/i, • ■ • ,Pmfm) , there exist polynomials hi, . . . , h v S 
P such that f = h%gi + ■ ■ ■ + h v g v and deg(ft.,) < mdo F (V) — 1 whenever 
higi ^ . 

(A2) For every non-zero polynomial f G ID J with a representation f — p\f\ + 
• • • + Pmfm and V = (pi/i, . . . ,p m fm) , there exist h\,...,h v e P such 
that f = h\g\ + • • • + h v g v and max{deg(/ij) | i G {1, . . . , v}, higi 7^ 0} = 
ind 0j p(7>)-l. 

Proof. First we show that (Ai) holds if G is an Of -border basis. The Subideal 
Border Division Algorithm computes a representation / = h\g\ + • • ■ + h v g v + 
Citif ai +- ■ ■+c ll t l j,f a with hi, ■ ■ ■ ,h v G P and c\, . . . , G K such that deg(/ij) < 
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indo F (P) — 1 for i = 1, . . . , v . Then c\t\f ai + • • • + c^t^fa^ = modulo /, and 
the hypothesis implies c\ = ■ ■ ■ = c M = . 

Next we prove that (Ai) implies {A 2 ) . If deg(/i;) < inde^T 3 ) — 1, then 
Prop. 13.41 2 shows that the Of -index of every representation of hi§i is at most 
deg(hi) + 1 and hence smaller than mdo F (V). By Prop. [3~4l 4. there has to be at 
least one number i G {1, . . . , v} such that deg(hi) — inde> F (V) — 1 . 

Finally, we assume (A2) and show the subideal border basis property. Let 
ci, . . . , c M € K satisfy citif ai H h c^t^fa^ G / n J. Then either / = Ci*i/ ai + 

• • • + c^t^fa equals the zero polynomial or not. In the latter case we apply (A2) 
and obtain a representation / = h\g\ + - ■ ■+h I/ g I/ with h%, . . . ,h v € P. Since / 7^ 0, 
we have max{deg(/ii) | i € {1, . . . , u}, higi 7^ 0} > 0. But m.do F ^P) — 1 = — 1 is 
in contradiction to the second part of (A2) . Hence we must have / = 0. Thus 
ID Jn (Of)k — 0, i.e. the set G is an 0j?-subideal border basis of /. □ 

Definition 3.7. In the setting of the algorithm, let Q = (gi, . . . ,<?„). Then the 
polynomial 

NRo F ,g(V) = ciii/ Ql H h c^t^f^ 

is called the normal remainder of the representation V — {pi fx, ■ ■ ■ ,p m fm) of / 
with respect to Q . 

Clearly, the normal remainder depends on the choice of the representation V . 
It has the following application. 

Corollary 3.8. In the setting of the algorithm, the residue classes of the elements 
of Of generate the image of the ideal J in P/(G) as a K -vector space. 

In other words, the residue class of every polynomial f € J can be represented 
as a K -linear combination of the residue classes {tif ai , . . . , i^fa } . Indeed, such 
a representation can be found by computing the normal remainder NRo Fi g(7- > ) for 
Q = (gi, ...,g v ) and the representation V = (pi/i, . . • ,p m fm) of f = pifi H h 

Pmfm • 

Proof. By the algorithm, every f € J can be represented in the form / — h\g\ + 

V Kg v + citif ai H h c^t^fa^ , where h\,...,h v € P and c 1( . . . , c M G K. 

Forming residue classes modulo (G) yields the claim. □ 

Our last corollary provides another motivation for the name "subideal border 
basis" . 

Corollary 3.9. In the setting of the algorithm, let G be an Of -subideal border 
basis of an ideal I C P . Then G generates the ideal I PI J . 

Proof. By definition, we have (g±, . . . , g v ) C Jn J. To prove the converse inclusion, 
let felDJ. Using the Subideal Border Division Algorithm, the polynomial / can 
be expanded as f — h 1 g 1 + - ■ ■+h u g v +c\t\f ai +- ■ ■+c lJi t^f a ^ , where hi, . . . ,h v G P 
and c%, . . . , Cfj, G K . This implies the equality of residue classes = / = c\iif ai + 

• • • + Cfjt^fa in P/I . By assumption, the residue classes i\f ai , ■ ■ ■ ji^fa form a 
if-vector space basis of (I + J)/ 1. Hence c\ = ■ ■ ■ = = 0, and the expansion 
of / yields / = h m + ■■■ + h v g v G (G) . □ 
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4. The Subideal Version of the BM- Algorithm 

Let K be a field, let P = K[x\, . . . , x„] be the polynomial ring in n indetermi- 
nates over K, equipped with the standard grading, and let T™ be the monoid 
of terms in P. Given a finite set of points X = {pi, . . . ,p s } C K n , we let 
eval : P — ► K s be the evaluation map eval(/) = (/(pi), ■ ■ ■ , f(Ps)) associated 
to X. It is easy to adjust the Buchberger-Moller Algorithm (BM- Algorithm) so 
that it computes a border basis of the vanishing ideal 

h = (feP\ f(pi) = ■■■ = f(p s ) = 0) = ker(eval) C P 

of X. Since we use a version which differs slightly from the standard formulation 
(see for instance [4] or [TO], Thm. 6.3.10), let us briefly recall its main steps. 

Algorithm 4.1. (BM- Algorithm for Border Bases) 

Let X = {pi, . . . ,p s } C K n be a set of points given by their coordinates, and let a 
be a degree compatible term ordering on T™ . The following instructions define an 
algorithm which computes the order ideal O a {T) = T n \ LT CT (ix) and the O a {Ix)- 
border basis G of Ix- 

Bl Let d= 0, O = {I}, G = 0, and M = (1, . . . , l) tr G Mat M (A). 

B2 Increase d by one and let L — [ti, . . . , tg\ be the list of all terms of degree d 

in dO , ordered decreasingly w.r.t. a . If L = 0, return (<D,G) and stop. 
B3 Form the matrix A — (eval(ti) | • • • | eval(tf) | M) and compute a matrix B 

whose rows are a basis of the kernel of A . 
B4 Reduce B to a matrix C — (cy ) € Matfc ^+ m (iT) in row echelon form. 
B5 For all j 6 {1, ...,£} such that there exists an i € {1, . . . , k} with pivot 

index v(i) = j , append the polynomial 

i e+m 

tj + c ij'tj' + Cij'Uji 

j'=j+l j'=t+l 

to the list G , where uy is the (j' — £) th element of O . 
B6 For all j = 1,1 — 1, ... , 1 such that the j column of C contains no pivot 
element, append the term tj as a new first element to O , append the column 
eval(ij) as a new first column to M. , and continue with step B2. 

The proof of this modified version is simply obtained by combining all the it- 
erations of the usual BM- Algorithm corresponding to terms of degree d into one 
"block" . The fact that we put the terms of degree d in dO into L in step B2 effects 
the computation of the entire border basis, rather than just the reduced tr-Grobncr 
basis of Ix (see [5], Thm. 3.3). A further elaboration is beyond the scope of the 
present paper and is left to the interested reader. 

Given X and a polynomial ideal J = (F) with F = {/i, . . . , f m } C P \ {0}, 
we know that the vanishing ideal Ix has an O a (I)f -subideal border basis. The 
following generalization of the BM-algorithm computes this subideal border basis. 

Algorithm 4.2. (Subideal Version of the BM- Algorithm) 

Let X = {p\, . . . ,ps\ C K n be a set of points given by their coordinates, let a be 
a degree compatible term ordering, and let F — {f%, . . . , f m } C P \ {0} be a set 
of polynomials which generate an ideal J = (F) . The following instructions define 
an algorithm which computes an F -order ideal 0„{I)f and the O a (I) f -subideal 
border basis G of Ix- 
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51 Let d = min{deg(/i),...,deg(/ m )} - 1, O f = 0, G = 0, and M e 
Mat s ,o(iT). 

52 Increase d by one. Let L = [tif ai , . . . ,tgf at \ be the list of all F -terms of 
degree d in F\JdOp , with their leading terms ordered decreasingly w.r.t. a . 
If then L — and d > max{deg(/i), . . . , deg(/ m )} , return (Of,G) and 
stop. 

53 Form the matrix A = (eval(fi/ Ql ) | ••• | eval(ttf ae ) \ M.) and compute a 
matrix B whose rows are a basis of the kernel of A . 

54 Reduce B to a matrix C = (c^) G Ma,tk.e+m(K) in reduced row echelon 
form. 

55 For all j G {1, ■■■,(} such that there exists an i € {1, . . . , k} with pivot 
index v(i) = j , append the polynomial 

i l+m 
tjfaij + c ij'tj'f aj , + Cij'Uj' 

j'=j+i f=e+i 

to the list G, where Uji is the (j' - £) th element of O f . 

56 For all j = £, I — 1, . . . , 1 such that the j column of C contains no pivot 
element, append the F -term tjf aj as a new first element to Of, append 
the column eval(tjf aj ) as a new first column to Ai , and continue with 
step S2. 

Proof. First we show finiteness. When a new degree is started in step S2, the 
matrix M. has m = H=Of columns where Of is the current list of -F-terms. In 
step S6 we enlarge M. by new first columns which are linearly independent of the 
other columns. This can happen only finitely many times. Eventually we arrive at a 
situation where all new columns eval(tj/ ai ) of A in step S3 are linearly dependent 
on the previous columns, and therefore the corresponding column of C contains a 
pivot element. Consequently no elements are appended to Of in that degree and 
we get L = in the next degree. Hence the algorithm stops. 

Now we show correctness. The columns of A are the evaluation vectors of F- 
terms whose leading terms are ordered decreasingly w.r.t. a . A row (en, . . . , cn +m ) 
of C corresponds to a linear combination of these F -terms whose evaluation vector 
is zero. Let gi,...,g^ be the polynomials given by these linear combinations of 
F-terms. Clearly, we have <?, S 1% n J . 

The evaluation vectors of the F-terms which are put into Of in step S6 are 
linearly independent of the evaluation vectors of the F -terms in the previous set Of 
since there is no linear relation leading to a pivot element in the corresponding 
column of C. Inductively it follows that the evaluation vectors of the F -terms 
in Of are always linearly independent. Henceforth the pivot elements of C are 
always in the "new" columns and the polynomials gi have degree d. By the way 
the algorithm proceeds, every F-term in the border of the final set Of appears in 
exactly one on the elements of G . All the other summands of a polynomial gi are 
in Of- Hence the final set G is an Of -subideal border prebasis. 

Furthermore, every F-term is either in Of or it is a multiple of an F-term 
in OOf (cf. Prop. [3T4l 5). In the latter case, its evaluation vector can be written 
as a linear combination of the evaluation vectors of the elements of Of ■ Thus 
the evaluation vectors of the elements of Of generate the space of all evaluation 
vectors of F-terms. Since they are linearly independent, they form a if -basis of 
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that space. Now we use the facts that evaluation yields an isomorphism of if -vector 
spaces eval : P/I — ► K s and that the residue classes of the F-terms generate the 
if -vector subspace (I + J)/I of P/I to conclude that the residue classes of the 
F -terms in the final set Of form a if -basis of (I + J)/I . □ 

Let us illustrate this algorithm by an example. 

Example 4.3. In the polynomial ring P = Q[x, y, z] , we consider the ideal J = (F) 
with F = {/i, f 2 } given by /1 = x 2 — 1 and f 2 = y — z. Let a = DegRevLex. 

We want to compute an Oi?-subideal border basis of the vanishing ideal of the 
point set X= {(1,1,1), (0,1,1), (1,1,0), (1,0,1)}. Notice that the first point of X 
lies on Z(fi,f 2 ), so that we should expect an F-order ideal consisting of three 
F -terms. Let us follow the steps of the algorithm. (We only list those steps in 
which something happens.) 

52 Let d = 1 and L = [y - z\. 

53 Form A = (0, 0, 1, -l) tr and compute B = (0) . (Thus C = B.) 
S6 Let O f = {y - z} and M = (0, 0, 1, -l) tr . 



S2 Let d = 2 and L = [x 2 — 1, x(y — z), y(y — z), z(y — z) 



S3 Compute A = 



-1 







1 

-1 







-1 





1 

-1/ 



and B = 



-1 



(Thus C=B. 

S5 The pivot indices v(l) = 2 and v{2) = 3 yield the set G — {31,32} with 
9i = x{y - z) - (y - z) and g 2 = y(y - z) + z(y - z) - (y - z) . 

/0 0\ 



S6 We obtain Of = {x 2 - 1, z(y - z),y - z) and M 



S2 Let d 



z), z 2 (y - z 



S3 Find A 



and B - 



We have L 


= [ 


x{x 


2 _ 


i), y{* 2 


-1) 


01- 
/o 














°\ 














-1 




















1 


\0 - 


-1 





-1 


-1 


-v 




/l 














°\ 




1 



















1 






















1 








-1 













1 













[p 








1 


-1 


0/ 



1 

\0 -1 

1), xz(y 




1 

-V 

-z), yz(y- 



(Thus C = B.) 



-- x(x 2 

97 = 



- 1), .94 = y(x 2 - 1), 

yz(y — z), and finally 



55 Here we obtain G = {gi, . . . ,g$} where 33 = 
g b = z(x 2 - 1), ,g 6 = xz(y - z) - z(y - z) 
38 = z 2 (y - z) - z(y - z). 

56 There are no new non-pivot indices. Hence O and M. are not changed. 
S2 We get L = and the algorithm stops. 

The result is the F-order ideal Of = {x 2 — 1, z(y — z), y — z} and the Of 
subideal border basis G = {gi, . . . , gs} of ix • 
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5. The Subideal Version of the AVI- Algorithm 



From here on we work in the polynomial ring P = R[xi, . . . , x n ] over the field of 
real numbers. We let X = {pi, ...,j),}c[-l,l]"C R™ be a finite set of points and 
e > t > two threshold numbers. (The number e can be thought of as a measure 
for error tolerance of the input data points X and r is used as a "minimum size" 
for acceptable leading coefficients of unitary polynomials.) 

Let us point out the following notational convention we are using: the "usual" 
norm of a polynomial / £ P is the Euclidean norm of its coefficient vector and is 
denoted by ||/||. By "unitary" we mean |j/|| = 1. In contrast, by ||/||i we mean 
the sum of the absolute values of the coefficients of / , and the term " || ||i -unitary" 
is to be interpreted accordingly. 

Furthermore, by eval : P — ► W we denote the evaluation map eval(/) = 
(/(pi), . . . ,f(p s )) associated to X. For the convenience of the reader, we briefly 
recall the basic structure of the Approximate Vanishing Ideal Algorithm (AVI- 
algorithm) from [5] . Notice that we skip several technical details and explicit error 
estimates. The goal of the AVI-algorithm is to compute an approximate border 
basis, a notion that is defined as follows. 

Definition 5.1. Let O — {tx, . . . , i M } C T" be an order ideal of terms, let dO — 
{bi, . . . , b v } be its border, and let G — {gi, . . . , g v } be an O-border prebasis of 
the ideal / = (<7i, . . . , gj) in P. Recall that this means that gj is of the form 

gj = bj - witn c v e R - 

For every pair such that bi, bj are neighbors in dO , we compute the normal 
remainder S[j = N~Ro,G{Sij) of the S-polynomial of gi and gj with respect to G . 
We say that G is an e -approximate border basis of the ideal I = (G) if we have 
||SV, || < e for all such pairs 

Moreover, the AVI-algorithm uses the concepts of approximate vanishing, ap- 
proximate kernel and stabilized reduced row echelon form, for which we refer to [5] , 
Sect. 2 and 3. 

Algorithm 5.2. (AVI- Algorithm) 

Let X = {pi, . . . ,p s } C [—1, 1]" C K™ be a set of points as above, and let a be a 
degree compatible term ordering. Consider the following sequence of instructions. 
Al Start with lists G = 0, O = [1], a matrix M = (1, . . . , l) tr G Mat s ,i(R), 
and d = . 

A2 Increase d by one and let L = [fi, . . . , tf\ be the list of all terms of degree d 
in dO , ordered decreasingly w.r.t. a . If L = , return the pair (0, G) and 
stop. 

A3 Form the matrix A = (eval(ti), . . . , eval(i^), M.) and calculate a matrix B 
whose rows are an ONB of the approximate kernel apker(yl, e) of A. 

A4 Compute the stabilized reduced row echelon form of B with respect to the 
given t . The result is a matrix C — (cy) € Matfc^_)_ m (R) such that cy = 
for j < . Here v(i) denotes the column index of the pivot element in 
the i row of C . 

A5 For all j € {1, . . . ,£} such that there exists an i € {1, . . . , k} with = j , 
append the polynomial 



I' 



t+m 




j'=j+l 



j'=t+l 
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to the list G , where Uj/ is the (J' ~ £) th element of O . 
A6 For all j — I, 1 — 1, . . . , 1 such that the j th column of C contains no pivot 

element, append the term tj as a new first element to O and append the 

column eval(ij) as a new first column to M. . 
A 7 Calculate a matrix B whose rows are an ONB of apker(.M, e) . 
A8 Repeat steps A4 - A7 until B is empty. Then continue with step A2. 

This is an algorithm which computes a pair (0, G) such that the following properties 
hold for the bounds 6 and r\ given in [5], Thm. 3.3. 

(a) The set G consists of unitary polynomials which vanish 5 -approximately at 
the points of X . 

(b) The set O — {t\, . . . , £ M } contains an order ideal of terms such that there 
is no unitary polynomial in (0}k which vanishes e -approximately on X. 

(c) The set G = {(1/ LC a (g)) g \ g € G} is an O -border prebasis. 

(d) The set G is an n -approximate border basis. 

Our main algorithm combines the techniques of this AVI-algorithm with the 
subideal version of the BM-algorithm presented above (see Alg. 14. 2|) . The result is 
an algorithm which computes an approximate subideal border basis. This notion 
is defined as follows. 

Definition 5.3. Let Of = {tif ai ,---,t^,f a } be an F-order ideal, let &Of = 
{bifPij-'-jbufp^} be its border, and let G — {gi, . . . ,g u } be an O^-subideal 
border prebasis. Recall that this means that c/j is of the form gj — bjffj j — 

J2i=i '•<..'</- wrt h Cij G R. 

For every pair such that bi, bj are neighbors in BOf , i.e. such that /Jj = j3j 
and bi,bj are neighbors in the usual sense, we compute the normal remainder 
S[j = NRojr.cC'Sy) of the S-polynomial of gi and gj with respect to G. We say 
that G is an e -approximate Of -subideal border basis if we have \\Sij\\ < e for all 
such pairs (i, j). 

Now we are ready to formulate and proof the main result of this section. 

Algorithm 5.4. (Subideal Version of the AVI- Algorithm) 

Let X = {pi,...,p s } C [—1,1]™ C R n be a set of points as above, let a be a 
degree compatible term ordering, and let F = {f\, . . . , f rn } C P \ {0} be a set of 
||i -unitary polynomials which generate an ideal J = (F) . Consider the following 
sequence of instructions. 

SA1 Let d = min{deg(/i),...,deg(/ m )} - I, O f = 0, G = 0, and M G 
Mat s ,oOFQ. 

SA2 Increase d by one. Let L = [t\f ai , . . . ,tif a( ] be the list of all F -terms of 
degree d in FUOOf , with their leading terms ordered decreasingly w.r.t. a . 
If then L = and d > max{deg(/i), . . . , deg(/ m )} , return (Of,G) and 
stop. 

SA3 Form the matrix A — (eval(ti/ ai ) | ••• | ev&\(tif ae ) | M) and compute a 
matrix B whose rows are an ONB of the approximate kernel of A . 

SA4 Compute the stabilized reduced row echelon form of B with respect to the 
given r. The result is a matrix C = (cij) € Matfc j £+ m (R) such that c\j — 
for j < v(i) . Here v{i) denotes the column index of the pivot element in 
the i th row of C . 
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SA5 For all j 6 {!,...,£} such that there exists an i € {1, . . . , k} with v(i) = j , 
append the polynomial 

l i+m 
tjfaj + Cij'tj'fctj, + Cij/Uji 

j'=j+l j'=t+l 

to the list G, where Uj> is the (j' - £) th element of O f . 
SA6 For all j = £, I — 1, . . . , 1 such that the j th column of C contains no pivot 

element, append the F -term tjf a . as a new first element to Of, append 

the column eval(tjf aj ) as a new first column to A4 . 
SA7 Calculate a matrix B whose rows are an ONB of apker(A^,e) . 
SA8 Repeat steps SA4 - SA7 until B is empty. Then continue with step A2. 

This is an algorithm which computes a pair (Of,G) with the following properties: 

(a) The set G consists of unitary polynomials which vanish S -approximately at 
the points of X. Here we can use S = s-s/p + tv(jj, + v)^fs. 

(b) The set Of contains an F -order ideal such that there is no unitary poly- 
nomial in (Of)k which vanishes e -approximately on X. 

(c) The set G = {(1/ LC CT (g)) g \ g € G} is an Of -subideal border prebasis. 

(d) The set G is an r\ -approximate subideal border basis for r] = 26 + 2vS 2 /j£ + 
2v5^/s/e. Here 7 denotes the smallest absolute value of the border F -term 
coefficient of one the polynomials gi . 

Proof. Large parts of this proof correspond exactly to the proof of the usual AVI- 
algorithm (see Thm. 3.2 in [5]). Therefore we will mainly point of the additional 
arguments necessary to show the subideal version. The finiteness proof is identical 
to the finiteness proof in the subideal version of the BM-algorithm 14.21 

For the proof of (a), we can proceed exactly as in the case of the usual AVI- 
algorithm. There is only one point where we have to provide a further argument: 
the norm of the evaluation vector of an i^-term is < y^s . To see this, we let tifj 
be an F-tevm and we write tifj = c^tk with Ck £ K and tk £ T" . Since fj 
is || ||i-unitary and X G [—1,1]" , we have || eval(ij/j)|| < Ylk \ Ck \ II eva l(^)l — 
ll/jlli\/s = Vs- 

Next we show (b). The columns of the final matrix M. are precisely the eval- 
uation vectors of the F -terms in Of- After the loop in steps SA4 - SA8, we 
have apker(Al) = {0}. Hence no unitary polynomial in (Of)k has an evalua- 
tion vector which is smaller than e. It remains to show that Of is an F -order 
ideal. Suppose that Ufj e Of and that XkUfj is put into Of- We have to prove 
that every .F-term t fj such that xgtfj = XkUfj is also contained in Of- In this 
case we have tj = X(t! and we want to show Xkt'fj € Of- For a contradiction, 
suppose that Xkt' fj is the border _F-term of some g € G. Since the evaluation 
vector of xgXkt' fj — x^tifj is not larger than ev&l(xkt' fj) , also this i^-term would 
be detected by the loop of steps SA4 - SA8 as the border F-term of an element 
of G. This contradicts XkXifj £ Of- 

To prove (c), it suffices to note that steps SA2 and SA5 make sure that the 
elements of G have the necessary form. Finally, claim (d) follows in exactly the 
same way as part (d) of [S], Thm. 3.3. □ 

Let us follow the steps of this algorithm in a concrete case which is a slightly 
perturbed version of Example 14.31 
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Example 5.5. In the ring P = M.[x,y,z) we consider the ideal J = (/i,/2) gen- 
erated by the || ||i -unitary polynomials f± = 0.5 y — 0.5 z and f 2 = 0.5 a; 2 — 0.5. 
Let a = DegRevLex, let e — 0.03, and let r = 0.001. We want to compute 
an approximate subideal border basis vanishing approximately at the points of 
X = {(1,1,1), (0,1,1), (1,1,0), (1,0,0.98), (0.98,0,1). 

Notice that the first point of X is contained in Z(hifz) and that the last two 
points of X differ by < e from one point (1, 0, 1) . Hence the approximate subideal 
border basis should correspond to three points outside Z(J) , and therefore we 
should expect to get an F -order ideal consisting of three F-terms. We follow the 
steps of the subideal version of the AVI-algorithm [5~4l 

SA2 Let d = 1 and L = [0.5 y - 0.5 z] . 

SA3 We compute A = (0, 0, 0.5, -0.49, -0.51) tr and B=(0). (Thus C = B .) 
SA6 Let O = {/J and M = (0, 0, 0.5, -0.49, -0.5) tr . 
SA2 Let d = 2 and L = [f 2 , xf u yh, zh] . 

( 


0.5 
-0.49 
-0.49 



SA3 We compute A 






0.5 









-0.4802 
-0.5 



and 



-0.5 



\-0.0198 

'0.0004 0.6755 -0.5089 -0.5068 -0.1667 
-0.3812 -0.3735 -0.3812 0.7548 
SA4 The stabilized reduced row echelon form of B is 

'0 0.3812 0.3735 0.3812 -0.7548 
y 0.5754 0.5811 -0.5754 / 

SA5 We get G = {gi,g 2 } with gi = 0.3812 xfi + 0.3735 yf x + 0.3812 zfi 
0.7548/i and g 2 = 0.5754 y/i + 0.5811 zfi - 0.5754/i . 



B = apker(^l, e) = 



C = 






0.5 
-0.49 
-0.5 / 

0.5068 
0.3812 



SA6 We find O = {/ 2 , z/i, fx] and M = 



( 
















-0.5 


-0.5 
























v- 


-0.02 





-0.02 







0.5 





-0.4802 
-0.49 





-0.5 

-0.5; 



z 2 fi] 









\ 


-0.5 














0.5 





-0.48 


-0.49 


-0.02 


-0.5 


-0.5 / 



SA2 Now let d = 3 and L = [xf 2 , yf 2 , zf 2 




SA3 A= 

-0.48 
-0.49 

and B, C are matrices of rank 6 which yield six further approximate subideal 

border basis elements. 
SA5 We obtain G = {gi,...,g s } with g 3 = xf 2 - 0.02 zfi, gi = 0.71 yf 2 - 

0.71/ 2 + 0.01 zh, 35 = 0.71 zj 2 - 0.71/ 2 , g e = 0.71 xzft - 0.7 zf u g 7 = 

yzh, and g 8 = 0.71 z 2 h - 0.7 zfi ■ 
SA5 Since there is no new non-pivot row index, Of and M. are not changed. 
SA2 In degree d — 4 we find L = and the algorithm stops. 

Hence the result is the F-order ideal O — {x 2 — 1, z(y — z), y — z} and the 
approximate Of -subideal border basis G = {g\, . . . , gg} . This confirms that there 
are three approximate zeros of G outside the two lines Z(h,f 2 )- 
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6. An Industrial Application 

In this section we apply the subideal version of the AVI-algorithm to an actual 
industrial problem which has been studied in the Algebraic Oil Research Project 
(see [2]). Viewed from a more general perspective, this application shows how one 
can carry out the suggestion made in the introduction, namely to use the subideal 
version of the AVI-algorithm to introduce knowledge about the nature of a physical 
system into the modeling process. 

Suppose that a multi-zone well consists of two zones A and B . During so-called 
commingled production, the two zones are interacting and influence each other. We 
have at our disposal time series of measured data such as pressures, temperatures, 
total production and valve positions. Moreover, during so-called test phases we can 
obtain time series of these data when only one of the two zones is producing. The 
following figure gives a schematic representation of the physical system and the 
measured variables. 
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Figure 1 . Schematic representation of a two-zone well 



The measured total production does not equal the sum of the individual pro- 
ductions calculated from the test data. The production allocation problem is to 
determine the contributions of the two zones to the total production when they are 
producing together. Here the contributions ca,cb of the zones are defined to be 
the part of the total production pab passing through the corresponding down-hole 
valves. Therefore we have pab = ca + c b , but there is no way of measuring ca 
and cb directly. In this sense the production allocation problem is to determine 
the contributions ca,cb from the measured data. 

Let the indeterminate xa represent the valve position of zone A and Xb the 
valve position of zone B . Here Xi = means that the valve is closed and xi = 
1 represents a fully opened valve position. Clearly, if valve A is closed, i.e. for 
points in the zero set Z((xa)) , there is no contribution from zone A, and likewise 
for B. By Hilbert's Nullstellensatz, this means that the polynomial pa modeling 
the production of zone A should be computed by using the subideal version of the 
AVI-algorithm with J = (xa) ■ Similarly, we want to force ps & ( x b) ■ 

Now we model the total production pab in the following way. We write pab = 
PA + Pb + Qab where qab is a polynomial which measures the interaction of the 
two zones. To compute qAB , we write it in the form 



Qab = !a ■ (xb ■ Pa) + Ib ■ (xa • Pb) 
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Notice that such a decomposition can be computed via the subideal version of the 
AVI-algorithm by applying it to the ideal J = (xbpa, xaPb) ■ The result will be a 
representation pab = PA + Pb + f ax BP A + IbXaPb ■ Here we observe that xa = 
implies pab — Pb because pa S (xa) ■ Analogously, we see that xb — implies 
Pab = PA , in accordance with the physical situation. 

The endresult of these computations is that the contributions of the two zones 
during commingled production can be computed from the equalities ca — (1 + 
/axb)pa and eg = (1 + ]bxa)Pb- At the same time we gain a detailed insight 
into the nature of the interactions by examining the structure of the polynomials 
fA, Jb ■ 
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